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LONG  TERM  GOALS 

Development  of  accurate  and  fast  advanced  statistical  and  dynamical  nonlinear  models  of  ocean 
surface  waves,  based  on  first  physical  principles,  which  will  improve  and  accelerate  both  long  term 
ocean  surface  waves  forecasts  and  prediction  of  strongly  coherent  events,  such  as  freak  waves,  tsunami 
and  wave-breaking. 

OBJECTIVES 

Finding  of  physically  correct  wind  input  term  for  Hasselmann  equation,  understanding  of  the  balance 
of  source  terms  in  Hasselmann  equation,  investigation  of  the  problem  of  interaction  of  different  scales 
on  the  ocean  surface  (sea  and  swell),  development  of  new  water  surface  analytic  models  and  methods 
of  their  solution. 

APPROACH 

Advanced  analytical  techniques:  Hamiltonian  formalism,  self-similar  solutions,  analytic  solutions  of 
integral  equations;  numerical  methods  for  solution  of  integral  and  pseudo-differential  equations; 
comparison  of  analytic  and  numerical  results  with  experimental  data 

WORK  COMPLETED 

•  Finding  of  the  new  wind  input  term  through  experimental,  theoretical  and  numerical 
approaches 

•  Theoretical  and  numerical  proof  of  nonlinear  interaction  term  domination  over  wind  input  and 
dissipation  terms  in  Hasselmann  equation. 

•  Detection  of  swell  feedback  by  sea  background  through  theory  and  experimental  data 

•  Design  of  new  method  for  numerical  integration  of  Hasselmann  equation 

•  Derivation  and  numerical  testing  of  one-dimensional  version  of  Zakharov  equation,  especially 
convenient  for  theoretical  and  numerical  study  of  freak  waves 
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RESULTS 


1.  New  wind  input  term  through  experimental,  theoretical  and  numerical  consideration 
Nowadays,  Hasselmann  equation  is  widely  accepted  model  of  wind-driven  seas  [1]: 


M  =  SN  +  ea±eN=S"+S"+s 

dt  dt  dk  dr 


(1) 


The  nonlinear  interaction  tenn  Snl  is  perfectly  known.  Unfortunately,  the  knowledge  of  wind  input  Sin 
and  white-capping  dissipation  Sdiss  tenns  is  poor.  Creation  of  reliable  and  well  justified  theory  of  Sin  is 
hindered  by  strong  turbulence  of  the  air  boundary  layer  over  the  sea  surface.  As  a  result,  there  is  a 
dozen  of  heuristic  models  of  Sjn  with  enonnous  scatter  of  parameters  [4]:  the  Donelan  model  [2] 

predicts  Sjn  about  five  times  higher  than  Hsiao-Shemdin  model  [3].  Situation  with  Sdiss  is  not  better - 
the  theory  is  not  developed,  the  experimental  data  are  far  from  being  complete.  The  forms  of  Sdjss  used 
in  operational  models  are  heuristic  and  badly  justified  [18]. 


We  managed  to  find  the  new  fonn  of  the  dissipation  tenn  Sdiss  through  experimental,  theoretical  and 
numerical  approaches.  Resio  and  Long,  2004  [5]  presented  data  of  wind  wave  spectra  from  several 
experimental  installations  in  non-traditional  fonn,  expressing  coefficient 


n  1  [if  2  Y/3  ~1/2 

/3  =  -a4luAcp)  -u0\g 

(2) 

in  Kolmogorov-Zakharov  energy  spectra 

F(k)  =  J3k~5'2 

(3) 

in  terms  of  specific  combination  of  velocities.  Here  a4  =  0.00553  ,  cp  -  spectral  peak  velocity, 
u0  =  1 .93  M/cex,  uA  -  wind  speed  at  the  height  equal  fixed  portion  A  =  0.065  of  the  spectral  peak 
wave  length  2n  /  kp  ,  where  kp  is  the  wave  number  of  the  spectral  peak.  It  is  remarkable  that 

practically  all  experimental  results  presented  in  [5]  exhibit  portion  of  the  spectrum  obeying  the  law  (3). 
Consider  Hasselmann  equation  (1)  in  "time-domain"  form 

3£(f’g)  =  s„  +  r(®,  0)E(m,  0)  (4) 

dt 

where  E(co,d)  =  N(k,t ) ,  k  =  °J—  is  the  energy  spectrum.  Using  dimensional  estimate 

g 


g 


(5) 


forcing  function  in  the  form  y(co,  9)  =  acol+s  f  (0)  and  self-similar  substitution 

E(o),0,t)  =  tp+qF((Otq) 

we  get  the  indices  of  self-similar  solution  (5): 

1  9q-l 

q  - -  p  =  — - 

5  +  1  2 

Comparison  of  the  indices  p  and  q  with  characteristics  of  the  regression  line  (2)  gives  the  index  of  the 
wind  source  5  =  4/3  and  parameterization  of  the  new  wind  input  term: 

Sin  oc  a 7/3  (6) 

Numerical  simulation  of  the  Hasselmann  equation  (4)  with  wind  input  term  (6),  presented  on  Fig.l  , 
reproduces  the  regression  line  (2)  with  high  accuracy  for  the  wind  speed  values  from  2  to  10  m/sec. 

We  have  found  new  wind  input  term  through  experimental,  theoretical  and  numerical  considerations, 
which  can  be  used  for  quality  prediction  improvement  by  operational  wave  forecasting  models. 


Fig.  1  Experimental,  theoretical  and  numerical  evidence,  presented  on  the  same  plot,  in  variables 


1 000//  as  a  function  of 


Dotted  line  —  experimental  regression  line  by  Resio  and  Long, 


2004  [5J;  dashed  line  —  theoretical  result,  reproducing  regression  line;  numerical  result:  crosses  — 
wind  speed  u=2.5  m/sec;  stars  —  wind  speed  u=5.0  m/sec;  rectangles  —  wind  speed  u=l 0.0  m/sec; 

triangles  —  wind  speed  u=20.0  m/sec. 
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2.  Energy  balance  in  wind-driven  seas  within  the  Hasselmann  equation 

The  problem  of  balance  of  different  terms  in  the  right-hand  side  of  the  Hasselmann  equation  (1)  is  a 
key  question  for  both  wind-wave  interaction  theory  and  modeling.  Today's  mainstream  emphasis  is  on 
developing  new  functions  Sin  and  Sdiss ,  but  not  correct  and  accurate  calculation  Snl .  Confusion  comes 
from  [6],  where  all  three  source  terms  in  (1)  have  been  compared  for  the  case  of  fully-developed 
(mature)  sea,  and  conclusion  has  been  made  that  terms  Sin  and  Sdiss  can  be  two-three  times  greater 

than  Snl .  We  show,  in  fact,  that  situation  is  opposite:  Snl  has  a  leading  role  in  balance  of  wind-driven 
seas.  The  analysis  is  based  on  decomposition  of  Snl  into  nonlinear  damping  TkNk  and  forcing  Fk 

Snl  =  Fk  -  TkNk  (7) 

where  I).-  positive  nonlinear  damping  decrement,  Nk—  spectral  density  of  wave  action.  Our  numerical 
and  analytical  results  show  that  TkNk  and  Fk  surpass  conventional  parameterizations  of  input  and 
dissipation  of  wind-driven  waves  by,  at  least,  one  order  of  magnitude,  see  Fig.2. 

An  additional  argumentation  is  presented  in  Fig.  3,  where  nonlinear  damping  decrement  TkNk  is 
compared  to  empirical  parameterizations  of  wind-wave  growth  given  by  different  authors  and  used  in 
the  today  wave  forecasting  models  as  an  option. 


x  10  7 


Fig.2  Decomposition  of  the  nonlinear  term  S,,/  (solid  line)  for  the  case  by  Komen  et  al.  (1984)  [6] 
into  nonlinear  forcing  (dashed)  and  damping  (dotted)  terms  (see  Eq.  7) 
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Fig. 3  Nonlinear  damping  coefficient  Tk  given  by  theoretical  estimate  and  by  the  numerical 

simulation  (dashed  and  solid  bold  curves,  correspondingly).  Conventional  dependencies  of  wind 
growth  increments  are  shown  by  thin  curves  with  authors'  names  in  the  legend 

Our  result  on  dominating  effect  of  nonlinear  interaction  on  wave  spectra  evolution  should  not  be 
interpreted  as  a  call  to  ignore  the  effects  of  wind  input  and  wave  dissipation.  The  leadership  of  Snl 
does  not  mean  that  we  disregard  wind  input  and  dissipation,  we  just  put  them  in  proper  place.  The 
strong  nonlinear  forcing  and  damping  that  compose  the  conservative  term  Snl  determine  strong 
relaxation  and  a  universality  of  spectral  shaping  due  to  inherent  wave  dynamics,  while  Sin  and  Sdiss  are 
responsible  for  growth  of  total  energy. 

Obtained  result  will  help  to  make  better  theoretical  estimates  of  solutions  of  Hasselmann  equation  and 
help  to  develop  simplified  approximations  to  source  terms. 

3.  Global  visual  observations  as  a  tool  for  discriminating  swell  and  wind  seas 

Today's  Voluntary  Observed  Ships  (VOS)  data  are  considered  as  a  source  of  information  on  climatic 
features  of  the  world  ocean  rather  than  experimental  background  of  studies  of  physical  mechanisms  of 
wind- wave  dynamics.  We  are  trying  to  change  this  tradition  by  applying  results  of  asymptotic  weakly 
turbulent  model  of  wind-driven  seas  [1],  [4],  [8]  to  extensive  data  base  of  the  Global  Atlas  of  Ocean 
Waves  [9],  [10]. 

The  global  visual  wave  observations  are  re-analyzed  within  the  theoretical  concept  of  self-similar 
wind-driven  seas  .  The  theoretical  criteria  of  discriminating  wind-driven  and  swell  seas  are  fonnulated 
and  shown  to  be  adequate  to  the  problem.  The  results  are  detailed  for  the  South  Pacifica,  which  wave 
climatology  based  on  VOS  data  is  well  studied  and  the  swell  component  is  well  pronounced.  The  core 
of  the  analysis  are  one-parametric  dependencies  "wave  height  -  wave  period"  Hs  =  CTZ  .  The 

reference  cases  [11],  [12], [13]  have  found,  correspondingly,  Z  =  5/3  ,  Z  =  3/2  and  Z  =  4/3  .  This  set 
of  exponents  Z  has  been  interpreted  recently  in  [1],  [4],  [8]  in  terms  of  spectral  fluxes  and  total  wave 
input.  An  alternative  reference  case  -  sea  swell  gives  an  opposite  signature  of  the  exponent  Z  =  -1  /  2 . 
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This  simple  criterion  was  used  and  appeared  to  be  robust  for  the  problem  of  swell-wind  sea 
discrimination. 

The  corresponding  exponent  Z  appears  to  be  slightly  higher  than  Z  =  -1  /  2 ,  that  implies  a  pumping  of 
sea  swell  by  wind-driven  sea  background  [17].  This  important  issue  is  considered  both  in  the  context 
of  methodology  of  obtaining  VOS  data  and  within  the  physics  of  the  mixed  sea.  This  result  contradicts 
to  commonly  accepted  vision  of  sea  swell  as  a  neutral  or  slightly  decaying  fraction  of  ocean  wave 
field. 

Prospects  of  further  study  are  quite  promising.  In  particular,  satellite  data  are  seen  to  be  used  for 
tracking  ocean  swell  and  for  studies  of  physical  mechanisms  of  its  evolution. 

4.  A  new  method  for  solving  the  exact  SnI  in  the  energy  balance  equation 


One  of  the  most  challenging  tasks  in  improving  the  wave  forecasting  models  is  solving  accurately  and 
in  a  reasonable  time  the  collision  integral  in  the  energy  balance  equation  (1).  Here  we  present  a  new 
method  for  solving  the  collision  integral  in  the  Hassehnann  equation  that  can  be  used  in  arbitrary 
depth.  The  method  has  in  principle  some  advantages  with  respect  to  previously  developed  methods 
such  as  Resio-Tracy  (RT)  [14],  Lavrenov  (LAV)  [15],  Masuda  (MAS)  [16]:  (i)  the  loci  are  straight  lines 
instead  of  eggs-shaped  in  RT  method;  (ii)  The  calculation  of  the  loci  is  exact  and  does  not  require 
iterative  procedure;  (iii)  it  does  not  require  interpolation  of  the  spectrum  in  the  frequency  domain,  but 
just  in  angle  (the  interpolation  is  therefore  ID  with  respect  to  2D  in  RT,  MAS  and  LA  V);  (iv)  The  MAS 
and  LA  V  methods  contain  some  singularities  that  have  to  be  treated  with  some  caution.  The  method 
presented  here  does  not  have  singularities  as  the  RT  method.  The  collision  integral  is 


cNx 

dt 


1  1 

-  +  - 


1 


N,  N2  N3 


N. 


4  J 


8(Ak)8(co)dk234 


The  difference  between  all  these  methods  is  the  way  the  8  -functions  in  frequency  are  removed.  The 
methods  RT,  LA  V  and  MAS  remove  the  momentum  8  -function  by  integration  over  wave  vector.  Then 
8  -function  needs  to  be  solved  numerically.  Here  we  follow  different  approach.  First  we  are  making 
transform  to  polar  coordinates  -  the  resulting  integral  is  a  multidimensional  one  -  three  integrals  in 
wave  numbers  (modulus)  and  three  in  angles.  Integration  over  two  angles  is  then  performed  using  the 
property  of  8  -unction  and  the  resulting  integral  is  transformed  in  angular  frequency  coordinates  (three 
integrals  in  frequency  and  one  in  angle).  After  using  the  property  of  8  -function  in  frequency,  the 
integration  is  removed  and  we  come  to  the  final  result: 


dNx 

dt 


1  1 

-  +  - 


N, 


1 

If 


'n4, 


k3co2co3co4 


©2©4 


sin(02  -  Of) 


dcoTdO, 


with  04  =  ®(k4  -  \kx  cos(#|)  +  k2  cos (02)  -  k3  cos(#3)|)  where  0  is  the  step-function.  A  similar  definition 
holds  for  ©2.  The  numerical  integration  of  SnI  consists  in  discretization  of  the  above  integral.  With 
respect  to  the  RT  method  the  loci  over  which  we  make  the  integration  are  straight  lines.  We  give  a 
number  of  graphical  examples  comparing  the  loci  of  integration  for  the  RT  method  and  the  method 
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here  developed.  In  Fig.4a  we  show  the  classical  egg-shaped  loci  obtained  from  the  RT  method.  In  Fig. 
4b  we  show  the  respective  loci  in  our  domain  of  integration 


-1.5  -1  -0.5  0  0.5  1  1.5 


1.5 


co4 


1 


0.5 


0 


1.5  1.6  1.7  1.8  1.9  2  2.1  2.2 


k4x 


co2 


Fig. 4.  a)  Egg-shaped  loci  in  the  RT  method,  b)  Loci  in  the  new  numerical  method 


We  developed  new  promising,  in  terms  of  accuracy  and  speed,  method  for  solving  Hsslmann  equation 
nonlinear  term. 


5.  New  canonical  equation  for  one-dimensional  surface  waves 
We  applied  canonical  transfonnation  to  the  water  wave  Hamiltonian 


H  = 


+  gqfdx  +  ^q{ig;  -{ky)fdx  +  ^q{ky){k{q{kv))  +  riii/xx}dx 


where 


\y-r] 

The  transfonnation  removes  not  only  cubic  nonlinear  tenns,  but  simplifies  drastically  fourth  order 
terms  in  the  Hamiltonian: 


H  =  jV (bkbdx  + — j|  b'\ 


-(bb'*  -b*b')-k\b\2 


dx. 


b  -  is  new  normal  canonical  variable.  This  transfonnation  explicitly  uses  the  fact  of  vanishing  exact 
four  waves  interaction  for  water  gravity  waves  for  2D  potential  fluid.  After  the  transfonnation  well- 
known  but  cumbersome  Zakharov  equation  is  drastically  simplified  and  can  be  written  in  X-space  in 
compact  way.  This  new  equation  is  very  suitable  for  both  analytic  study  and  numerical  simulation: 


.db 

i  — 

dt 


d>kb  + 


d 

dx 


b-k(\b'\2)-^(b'k(\b\2)) 

OX 
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Also,  this  equation  is  convinient  for  experimentalists.  The  simplest  solution  of  the  equation  is  the 


monochromatic  wave  b(x )  =  B0e 


i(kQx-w0t) 


1  , 

with  the  frequency  co0  =  cok  +  —  k0\  B0  \~ .  Growth-rate  of 


modulation  instability  is  also  calculated: 


r\  =T“Gr(1“6//2^2 

8  k0 


/r(k0-^)2-  — 


In  the  framework  of  this  “improved”  Zakharov  equation  we  have  perfonned  numerical  simulation  of 
freak- wave  formation  from  the  initialy  unifonn  water  waves,  see  Fig.  I 


Re(b(x,t))  and  |b(x,t| 


Fig.  5  b(x)  and  Re{b{x)). 

We  developed  simple  new  equation  for  ID  nonlinear  surface  waves,  which  is  convenient  for  both 
theoretical  study  and  numerical  simulation  of  wave  propagation  in  narrow  experimental  tanks. 

IMPACT/APPLICATIONS 

•  Obtained  results  are  expected  to  have  significant  impact  on  operational  wave  models  through 
improved  wind  input  term  and,  as  a  consequence,  better  prediction  of  ocean  surface  wave 
evolution. 

•  Correct  understanding  of  Hassehnann  equation  source  terms  balance  will  help  to  develop  better 
theoretical  estimates  and  approximations  to  the  source  terms. 

•  Understanding  of  swell  -  wind  sea  interaction  has  to  help  to  develop  approximations  to 
nonlinear  source  term  and,  as  a  result,  build  simplified  theoretical  and  fast  numerical  models  of 
nonlinear  interaction  of  surface  waves. 

•  New  method  of  numerical  integration  of  Hassehnann  equation  will  help  to  develop  faster 
working  operational  wave  prediction  models. 
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•  New  one-dimensional  nonlinear  equation  for  surface  waves  will  help  to  carry  out  more 

effective  theoretical  analysis  for  freak  waves  and  faster  numerical  simulation  of  the  processes 
of  wave  propagation  in  narrow  experimental  tanks. 
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